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Abstract In this work we show results for light nuclear systems and small clusters of helium atoms 
using the hyperspherical harmonics basis. We use the basis without previous symmetrization or an- 
tisymmetrization of the state. After the diagonalization of the Hamiltonian matrix, the eigenvectors 
have well defined symmetry under particle permutation and the identification of the physical states 
is possible. We show results for systems composed up to six particles. As an example of a fermionic 
system, we consider a nucleon system interacting through the Volkov potential, used many times in 
the literature. For the case of bosons, we consider helium atoms interacting through a potential model 
which does not present a strong repulsion at short distances. We have used an attractive gaussian po- 
tential to reproduce the values of the dimer binding energy, the atom-atom scattering length, and the 
effective range obtained with one of the most widely used He-He interaction, the LM2M2 potential. In 
addition, we include a repulsive hypercentral three-body force to reproduce the trimer binding energy. 

1 Introduction 

The Harmonic Hyperspherical (HH) method is extensively used in the description of few-body systems. 
For example the HH method has been applied to describe bound states of A = 3, 4 nuclei (for a recent 
review see Refs. [2I2])- In these applications the HH basis elements, extended to spin and isospin 
degrees of freedom, have been combined in order to construct antisymmetric basis functions; in fact, 
the HH functions, as normally defined, do not have well defined properties under particle permutation, 
but several schemes have been proposed to construct HH functions with an arbitrary permutational 
symmetry, see Refs. [3UHEH]. 
. All of the proposed symmetrization schemes share an increasing computational difficulty as the 

number of particles A increases; to cope with this issue, the authors proposed in Ref. [7] to renounce to 
the symmetrization step. If the Hamiltonian commutes with the group of permutations of A objects, 
Sa, the eigenvectors can be organized in accordance with the irreducible representations of Sa', in fact, 
if there is no more degenerancy, the eigenvectors have a well defined permutation symmetry. After the 
identification of the eigenvectors belonging to the desired symmetry, the corresponding energies are 
variational estimates. The disadvantage of this method results in the large dimension of the matrices 

Presented at the Sixth Workshop on the Critical Stability of Quantum Few-Body Systems, Erice, Sicily, October 
2011 

M. Gattobigio 

Universite de Nice-Sophia Antipolis, Institut Non-Lineaire de Nice, CNRS, 1361 route des Lucioles, 06560 
Valbonne, France 

E-mail: mario.gattobigio@inln.cnrs.fr 
A. Kievsky - M. Viviani 

Istituto Nazionale di Fisica Nucleare, Largo Pontecorvo 2, 56127 Pisa, Italy 



2 



to be diagonalized. However, at present, different techniques are available to treat (at least partially) 
this problem. 

In order to show the main characteristics of this method, we will discuss results for bound states 
up to six particles interacting through a central potential in two different systems: (i) a nucleon system 
interacting via the Volkov potential, used many times in the literature [51I51[51 HU1[TT1[T^] . and thus useful 
to test our approach, and (ii) a systems composed by helium atoms interacting through a soft-core 
potential. The ab-initio-helium potentials have a strong repulsion at small distances which makes 
calculations quite difficult; few calculations exist on clusters of helium with these potentials [I51I14IIT5] . 
On the other hand, descriptions of few-atoms systems using soft-core potentials are currently operated 
(see for example Ref. HMIZ1II1])- 

The paper is organized as follows. In Sect. [5] a brief description of the method is given. In Sectf31 
applications of the method to a system of nucleons and helium atoms are shown. The conclusions are 
given in the Sect. |U 



2 The unsymmetrized HH expansion 

In the present section we give a brief description of the HH basis showing some properties of the basis 
that allow to use unsymmetrized basis elements to describe a system of identical particles. 



2.1 The HH basis set 

Following Refs. [lllH8U 7]. we start with the definition of the Jacobi coordinates for an equal mass A 
body system, with Cartesian coordinates ri . . . ta 



x^_ j+1 = ^^ T (r J - +1 -X i ), j = l,...,N. (1) 

with Xj = Si=i r i/j- For a gi ven se t of Jacobi coordinates xi,...,xjv, we can introduce the hy- 
perspherical coordinates. A useful tool to represent hyperspherical coordinates is the hyperspherical 
tree. This is a rooted-binary tree whose leaves represent the modules of Jacobi coordinates. Once we 
introduce the hyperradius, 

/ N \ 1/2 / A , 1/2 / „ A \ 1/2 

p=(E«?) = 2 5>~ x)2 =(i£(^) 2 ) , (2) 

^ i=l ' ^ i=l ' ^ j>i ' 

the modulus of the Jacobi coordinates live in a (N — l)-sphere of radius p and we can introduce 
N — 1 hyperangles to express the Jacobi coordinates as a function of the hyperradius. The choice is 
not unique, and different choices are represented by different hyperspherical trees [191 20 . The relation 
between a tree and the corresponding Jacobi coordinates is the following; for each tree's node, labelled 
by a, we have an hyperangle <f> a . The rule to reconstruct the value of a Jacobi coordinate modululs 
reads: start from the root node, and look for the path leading to the leaf corresponding to the Jacobi 
coordinate; for each branch turning toward the left (right) we multiply the hyperradius p for cosine 
(sine) of the hyperangle attached to the branching point. As an example, in Eq. ([3]) we have a tree 
choice for A — 5 with the corresponding relations between Jacobi and hyperspherical coordinates 

£4 £3 aS2 x\ 

X\ = p sin <^>4 sin <p2 
X2 = p sin <j)i cos <p2 
X3 = p cos (j)4 sin 03 

X4 = p COS 04 COS 03 , 




(3) 



K 



K 



Fig. 1 In the left panel we have drawn the standard hyperspherical tree; this is the one used in the standard 
definition of the basis, and the one used to calculate the two-body potential between particles at ri and r-2. In 
the right panel we have drawn the non-standard tree, used to calculate the three-body force between particles 
at n, r2, and r3. 



Different trees have different topologies; given a node a, the left (right) branch connects the node to 

a sub-binary tree made up of 
useful topological numbers as 



a sub-binary tree made up of N„ nodes and L 1 ^ leaves. We can use this information to construct 



and 



Sa=K+ l -L[ l+ l - 



(4) 



(5) 



The set of the hyperangles together with the direction of the Jacobi coordinates x, = [tpi, 6j) form 
the hyperangular coordinates 

Q N = (xi,...,Xjv, <Aa, •••,<£#•)• (6) 

in terms of which the HH functions D^fR:] (^jv) are defined. The subscript [K] stands for the set of 
(3JV — l)-quantum numbers l±, . . . , In, mi, . . . , mjsr, K2, ■ ■ ■ , -Kjv, with Kx = K the grand-angular 
momentum. They can be expressed in terms of the usual harmonic functions Yj m (x) and of the Jacobi 
polynomials P£' b (z) 



N 



j=i 



LM 



n "2 



.atE nodes 



with 



V 



f>a) =Mn 



■(cos0 o ) A -(sin0 a ) A «P, 



(COS 20 Q ) , 



where we have defined 



The normalization factor reads 



Af a0 



'2(2n + a + /3 + l)n! T(n + a + /3 + 1) 
r(n + a + l)r(n + /3 + l) 



(7) 

(8) 
(9) 

(10) 



With the above definitions, the HH functions have well defined total orbital angular momentum L and 
z-projection M. The standard choice of hyperspherical coordinates, and of the corresponding HH, is 
represented in the left panel of Fig. Q] this is the one we use as our basis set. 
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2.2 Rotation matrices between HH basis elements of different Jacobi coordinates 

Here we are interested in a particular set of coefficients relating the reference HH basis to a basis 
in which the ordering of two adjacent particles have been transposed. In the transposition between 
particles j,j + 1, only the Jacobi vectors Xj and Xj+i, with i — N — j + 1, are different. We label them 
x£ and x' i+1 , and explicitly they are 



i + iy 2(j + 1) 

Xj H ; X i+ i 

J J 



, y/U + l) 2 - 2(j + 1) 1 



(11) 



with i — 1, . . . , N — 1. The corresponding moduli verify x'j + x 1 i+1 — xf + xf +1 . Let us call yffi 
the HH basis element constructed in terms of a set of Jacobi coordinates in which the z-th and i + 1-th 
Jacobi vectors are given from Eq. ( HID with all the other vectors equal to the original ones (transposed 
basis). The coefficients 

ti,LM _ i jd [-y,LM / r> \T*^jLM/ 



^miK*] = J dQ N [y^(Q N )}*y^^ N ) , (12) 

are the matrix elements of a matrix Af M that allows to express the transposed HH basis elements 
in terms of the reference basis. The total angular momentum as well as the grand angular quantum 
number K are conserved in the above integral (K = K'). The coefficients -A^k^^ form a very sparse 
matrix and they can be calculated analytically using angular and T-coupling coefficients (Kil'dyushov 
coefficients) and the Raynal-Revai matrix elements [TTll21[[l~S] . 

We are now interested in obtaining the rotation coefficients between the reference HH basis and a 
basis in which the last Jacobi vector is defined as x^, = r, — , without loosing generality we consider 
j > i. A generic rotation coefficient of this kind can be constructed as successive products of the 
AicHK'] coefficients. Defining ' (^%) the HH basis element constructed in terms of a set of Jacobi 
coordinates in which the iV-th Jacobi vector is defined x^ = r 3 — r<, the rotation coefficient relating 
this basis to the reference basis can be given in the following form 

- J M&ffimryffiifft) = [^ M --^t M ] m • ( l3 ) 

The particular values of the indices ii, . . . , i n , labelling the matrices Af ± M , . . . , , depend on the 
pair The number of factors cannot be greater than 2(j — 2) and it increases, at maximum, by 

two units from j to j + 1 . The matrix 

Bg* = j*™-.j*i M , (14) 

is written as a product of the sparse matrices A\ M 's, a property which is particularly well suited for 
a numerical implementation of the potential energy matrix. 



2.3 The two-body and three-body potential energy matrices 

We consider the potential energy of an A-body system constructed in terms of two-body interactions 

V = J2V(i,j) . (15) 

i<j 

In the case of a central two-body interaction, its matrix elements in terms of the HH basis are 

V [K][KI] ( P ) = E(^Y(^)I^(^J)I^(^)) ■ (16) 

i<j 
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In each element (y(jq\V(i,j)\yhM) the integral is understood on all the hyperangular variables and 
depends parametrically on p. Explicitly, for the pair (1,2), the matrix elements of the matrix Vi2(p) 
are 

VmUip) = (y$(a N )\v(i,2)\y™(f2 N )) = 

Si u i[ ■ ■ ■ 5i n ,v n 5l 2 ,l' 2 ■ ■ ■ 5l n ,l' n 5k 2 ,k' 2 • • • &k n ,k' n (17) 
x / ^(cos^sin^) 2 ^V^-'iMVipcosM (V/*- 1 ^). 



(18) 



Using the rotation coefficients, a general term of the potential V(i,j) results 



V, 



(hi) 



[K][K 



or, in matrix notation, 



LM j^ij,LM 



[ x ] - [ ^] [ ^]<^](^)l^(i,i)l^](^)) 



[K»][K>»] 



The complete potential matrix energy results 



(19) 



(20) 



i<j 



Each term of the sum in Eq. (|20[) results in a product of sparse matrices, a property which allows 
an efficient implementation of matrix- vector product, key ingredient in the solution of the Schrodinger 
equation using iterative methods. 

The three-body force used in the present work depends on the hyperradius pijk of a triplet of 
particles r^,r,-,rfc. For an A-body systems, there are (~5) three-body terms 



V® = W( Pijk ), 

i<j<k 



(21) 



and one of them is the force between the triplet ri, r2, r3 for which we have p\ 2 3 = x 2 N + x 2 N _ 1 . This 
term can be easily calculated on a hyperspherical-basis set relative to an non-standard hyperspherical 
tree with the branches attached to leaves xn and xn-i going to the same node, as shown in the 
right panel of Fig. [T] The transition between this tree and the standard tree is simply given by the 
T-coefficients 



Xjv,ijy In-1 Kjv- 



Xv.iju 



X,V-l,iiV-l K N-2 




njy — l ^jv — l K 

n»-l 



K, n,v 




(22) 



K, njv 



or 



y [K ] wn) = r ««-i k y { m ™ 



(23) 



where all the variable with the tilde refer to the non-standard tree. In fact, with this choice we simply 
have 

P123 = pcoscj) N , (24) 
and the fixed-rho matrix elements of the matrix W\23(p) are 
yLMrn \\ixrr„ \hiiM / 



(y^ ] (n N )\w( Pl23 )\y^(n N )) = 

fy.ii • • • 5i> NtlN S L > 2iL2 ■ ■ ■ S l ^ l 5 M ',m 5 k 2 ,k 2 ■ ■ ■ s k> n _ 1 ,k n - 1 : 



(25) 



(cos<p N f K (sm<j) N ) SK d^ N r K ; 



A 



(<f> N )W(p cos 0at) , 
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where Ck and Sk are the topological quantum numbers relative to the grand-angular-i^ root node. 
In practice the matrix is extremely sparse, and it is diagonal on all quantum numbers but the grand- 
angular momentum. 

The three-body force matrix in the standard basis is obtained by means of the 7~-coefficients 

{y^ ] {fiN)\w{p 12i )\y^{Q N )) = 

E n -a I<N _ 2 a lN _ 1 a lN a KN _ 2 a lN _ iaiN , LM , ^ n| H/ , s ]v LM f fj \\ (26) 

Sjv-i 

which for all practical purposes reduces to a product of sparse matrices. 

In order to calculate the other terms of the three-body force, we use the matrices A^ M , defined in 
Eq. (fT2|) . that transpose particles; with a suitable product of these sparse matrices 

V™ = A$™---A™, (27) 

we can permute the particles in such a way that xjv = i - ,; — r^, and xjv_i = 2/v / 3(i"fc — (r^ + Yj)/2), 
and p?- fc = x 2 N _ 1 + x 2 N , and the total three-body force reads 

V^= £ [D^YWiM^f- (28) 

Kj<k 



3 Applications of the HH expansion up to six particles 



In this section we present results for A — 3 — 6 systems obtained by a direct diagonalization of the 
Hamiltonian of the system. The corresponding Hamiltonian matrix is obtained using the following 
orthonormal basis 



(pf2\m[K}} = [p 



- ( «(«+l)/2. 



(a + m)\ 



(29) 



where Lm\f3p) is a Laguerre polynomial with a = 3N — 1 and /3 a variational non-linear parameter. 
The matrix elements of the Hamiltonian are obtained after integrations in the p, f2 spaces. They depend 
on the indices m, m' and [K], [K'] as follows 



(m' [K']\H\m[K]) 



- K(K + 3N- 2)T^J6 [K , ][K] 



+E 
E 

i<j<k 



^ij,LM *oij,LM -frin^ir. 

rin^iirnntirn " twin 



E B \k][K"]"[K'"][K'] r [K"][K'"\ 
[K"][K'"\ 

U [K\[K"\ U [K>"][K>] W [K"][K>"\ 

\K"]\K">] 



(30) 



The matrices T^ 1 - 1 and have an analytical form and are given in Ref. [TT]. The matrix ele- 
ments are obtained after integrating the matrix V\2 (p) in p-space whereas the matrix ele- 
ments W^'™ K ,^ are obtained after integration of the matrix 1^123 (p) (we will call the corresponding 
matrices V12 and W123, respectively). Introducing the diagonal matrix D such that ([K 1 ] \D\[K]) = 
8\k] ,[K']K{K +3N— 2), and the identity matrix I in i^-space, we can rewrite the Hamiltonian schemat- 
ically as 

H = -^-{ {1) T ® / + (2) T ® D) + ^[BfFY V 12 B% M + ]* ^123 V%f , (31) 

ij ijk 

in which the tensor product character of the kinetic energy is explicitly given. 
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Fig. 2 In the left panel we have the convergence of the 3 H and 3 He binding energies as a function of K. The 
excited state of the alpha particle 4 He* is also shown. In the right panel we have the convergence of the 4 He, 
6 He and 6 Li binding energies as a function of K. 



3.1 nuclear system 

As a first example we consider e nuclear system interaction through a simple two-body potential, the 
Volkov potential 

V{r) = V R e- r2 ' R * + V A e- r2 l R - , (32) 

with Vr = 144.86 MeV, Ri = 0.82 fm, V A = -83.34 MeV, and R 2 = 1.6 fm. The nucleons are 
considered to have the same mass chosen to be equal to the reference mass m and corresponding to 
h 2 /m = 41.47 MeVfm 2 . With this parametrization of the potential, the two-nucleon system has a 
binding energy E 2 n — 0.54592 MeV and a scattering length a 2 N — 10.082 fm. This potential has 
been used several times in the literature making its use very useful to compare different methods [SJ 
8 , 9 , lU] . The use of central potentials in general produces too much binding, in particular the A = 5 
system results bounded. Conversely, the use of the s-wave version of the potential produces a spectrum 
much closer to the experimental situation. This is a direct consequence of the weakness of the nuclear 
interaction in p-waves. Accordingly, we analyze this version of the potential, the s-wave projected 
potential. The results are obtained after a direct diagonalization of the Hamiltonian matrix of Eq. (1301) 
including m max + 1 Laguerre polynomials with a fix value of /3, and all HH states corresponding to 
maximum value of the grand angular momentum K max . The scale parameter /3 can be used as a non- 
linear parameter to study the convergence in the index m = 0, 1, . . . 

5 T^max j with Tn rnax the maximum 

value considered. We found that 20 Laguerre polynomials (with proper values of /3) were sufficient for 
an accuracy of 0.1% in the calculated eigenvalues. 

The results of the present analysis are given in Fig. [5] where the convergence of the A = 3 — 6 
binding energies are given as a function of K. In the left panel of Fig. [5]the convergence for the excited 
state 4 He* of the a particle is also shown. For A — 3, 4 a very extended HH expansion has been used 
with the maximum value of K — 80 and K — 40 respectively. For A = 3, the obtained results are 
8.431 MeV and 7.725 MeV for 3 H and 3 He respectively. For A — 4, the ground state binding energy 
converges at the level of 1 — 2 keV for K max = 40. The convergence of the excited state 4 He* has been 
estimated at the level of 50 keV. Though the convergence was not completely achieved, the description 
is close to the experimental observation of a + resonance between the two thresholds and centered 
395 keV above the p- 3 H threshold. Besides its simplicity, the s-wave potential describes the A — 3,4 
system in reasonable agreement with the experiment. 

For the A = 6 system a maximum value of K — 22 has been used which greatly improve previous 
attemps in using the HH basis in A = 6 systems @BB]. The obtained results are 33.016 MeV and 
32.087 MeV for 6 He and 6 Li respectively. It should be noticed that these states belong to the mixed 
symmetries [4 2] (without the Coulomb interaction). When the Coulomb interaction between two 
nucleons is included the state belong to the symmetry [2] <g) [2 2 ] and when it is included between three 
nucleons the state belongs to the symmetry [2 1] ® [21]. These states are embedded in a very dense 
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spectrum. In order to follow these state in the projected Lanczos method a projection-purification 
procedure is performed. 

3.2 atomic system 

As an example of an atomic systems we describe a system of 4 He atoms up to six atoms. The 4 He- 4 He 
interaction presents a strong repulsion at short distances, below 5 a.u. This characteristic makes it 
difficult a detailed description of the system with more than four atoms. Accordingly, we study small 
clusters of helium interacting through a soft-core two- and three-body potentials. Following Refs. [2"2l 
IT71H8] we use the gaussian two-body potential 

V(r) = V e~ r2 / R2 , (33) 

with V = -1.227 K and R = 10.03 a.u.. In the following we use h 2 /m = 43.281307 (a.u.) 2 K. This 
parametrization of the two-body potential approximately reproduces the dimer binding energy , the 
atom-atom scattering length ao and the effective range ro given by the LM2M2 potential. Specifically, 
the results for the gaussian potential are E2 — —1.296 mK, ao = 189.95 a.u. and r$ = 13.85 a.u., to be 
compared to the corresponding LM2M2 values E2 = —1.302 mK, ao = 189.05 a.u. and ro = 13.84 a.u.. 
As shown in Ref. [T7], the use of the gaussian potential in the three-atom system produces a ground 
state binding energy E^ — 150.4 mK, which is appreciable bigger than the LM2M2 helium trimer 
ground state binding energy of 126.4 mK. In order to have a closer description to the ^4 = 3 system 
obtained with the LM2M2 potential, we introduce the following three-body interaction 

W(p ijk ) = Wo e" 2 "WPo t (34) 

where p^ fe = |(ry + r 2 k + rj^) is the three-body hyperradius in terms of the distances of the three 
interacting particles. Moreover, the strength Wo is fixed to reproduce the LM2M2 helium trimer binding 
energy of 126.4 mK. In Ref. [TH] a detailed analysis of this force has been performed by varying the 
range po between 4 and 16 a.u.. Here we present results for small clusters, up to A = 6, formed by 
atoms of 4 He using the soft two-body force plus the hyperradial three-body force with parameters 
W ,po = 0.422K,14.0 a.u.. 

The results are collected in Figs. 13141 where we show the convergence in terms of K of the ground 
state and first excited state of the bosonic helium clusters. Starting from A — 3 the bosonic spectrum 
is formed by two states, one deep and one shallow close to the threshold formed by the A — 1 system 
with one atom far away. The calculations have been performed up to K = 40 in A — 4, K = 24 in 
A = 5 and K = 22 in A = 6. From the figure we can observe that the ground state binding energy, 
E { °\ has a very fast convergence in terms of K. The convergence of the Ev' is much slower than for 
the ground state, however with the extended based used it has been determined with an accuracy well 
below 1%. The results confirm previous analyses in the four body sector that the lower Efimov state in 
the A = 3 system produces two bound states, one deep and one shallow. Here, we have extended this 
observation up to the A = 6 system. Specifically we have obtained the following ground state energies: 
e[ 0) = 568.8 K, E { 5 0) = 1326.6 K,and £^ 0) = 2338.9 K, and first excited state energies: Ef ] = 129.0 

K, EP = 574.9 K,and E^ 0) = 1351.6 K. It is interesting to compare the results obtained using the 
soft-core representation of the LM2M2 potential with the results of Refs. [TSlITi] obtained using the 
original LM2M2 interaction. For the ground state the agreement is around 2% for A = 4, 5 and around 
1% for A — 6. The agreement is worst for the excited state, however the results from Ref. are 
obtained using approximate solutions of the adiabatic hyperspherical equations. The recent, and very 
accurate, results of Ref. [TS] for A — 4 (£4°' = 558.98 K and E^ — 127.33 K) shows a good agreement 
in particular for the first excited state. 

Finally it is possible to analyze the ratios Ej^ j E^-x- In the case of Efimov physics these ratios 
present and universal character. The He-He potential it is not located exactly at the unitary limit 
(infinite value of ao) but it is close to it. Using the soft potential models these ratios are: E± /E^ = 
1.020, E^/Ef 1 = 1.011 and E^/Ef 1 = 1.018. The ratios between the trimer ground state and the 

ground states of the A = 4,5,6 systems are Ef ] /Ef ] = 4.5, E { ° ] / Ef ] = 10.5 and E^/E^ = 18.5, 
respectively. These ratios are in good agreement with those given in Refs. [2"3"ll2~4"ll2"5"] . 
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The overall agreement of the A = 4, 5, 6 systems between LM2M2 and the soft potential model 
gives a further indication that these systems are a nice realization of which it is called Efimov physics. 
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Fig. 3 The trimer bound state and tetramer first excited state (left panel) and tetramer bound state and 
pentamer first excited state (right panel) as a functions of K. 
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Fig. 4 The pentamer bound state and hexamer first excited state (left pannel) and hexamer bound state (right 
pannel) as a functions of K. 



4 Conclusions 

In this work we have shown results using the HH expansion in the description of a A-body system with 
A = 3 — 6. The basis has not been symmetrized or antisymmetrized as required by a system of identical 
particles. However, the eigenvectors of the Hamiltonian have well defined permutation symmetry. The 
benefit of the direct use of the HH basis is based on a particular simple form used to represent the 
potential energy. We have limited the analysis to consider a central potential. In a first example we 
have describe a system of nucleons interacting through the Volkov potential, used several times in 
the literature. Though the use of a central potential leads to an unrealistic description of the light 
nuclei structure, the study has served to analyze the characteristic of the method: the capability of the 
diagonalization procedure to construct the proper symmetry of the state and the particular structure, 
in terms of products of sparse matrices, of the Hamiltonian matrix. The success of this study makes 
feasible the extension of the method to treat interactions depending on spin and isospin degrees of 
freedom as the realistic NN potentials. A preliminary analysis in this direction has been done 26 . 

In a second example we have studied the possibility of calculating bound and excited states in a 
bosonic system consisting of helium atoms interacting through soft two- and three-body forces. The 
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potential model has been adjusted to approximate the description of small helium clusters interacting 
through one of the realistic helium-helium interactions, the LM2M2 potential. After the direct diag- 
onalization of the A-body system we have observed that clusters with A = 3,4,5,6 atoms present a 
deep bound state and a shallow bound state just below the energy of the A — 1 system. Since the He-He 
potential predicts a large two-body scattering length we have studied the universal ratios Ej^/Ef^_\y 
These ratios have been studied in detail in the A — 4 case (see Refs. ;271124|). Estimates have been 
obtained also for bigger systems |25) . Our calculations, obtained for one particular value of the ratio 
ro/a, are in agreement with those ones. An analysis of the universal ratios as a — > oo is at present 
under way. 
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